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A novel approach for the refinement of finite-element-based analytical models of flexible structures is presented. 
The proposed approach models the possible refinements in the mass, damping, and stiffness matrices of the finite 
element model in the form of a constant gain feedback with acceleration, velocity, and displacement measurements, 
respectively. Once the free elements of the structural matrices have been defined, the problem of model refinement 
reduces to obtaining position, velocity, and acceleration gain matrices with appropriate sparsity that reassign 
a desired subset of the eigenvalues of the model, along with partial mode shapes, from their baseline values to 
t ose obtained from system identification test data. A sequential procedure is used to assign one conjugate pair 
of eigenvalues at each step using symmetric output feedback gain matrices, and the eigenvectors are partially 
assigned, while ensuring that the eigenvalues assigned in the previous steps are not disturbed. The procedure 
can also impose that gain matrices be dissipative to guarantee the stability of the refined model. A numerical 
example, involving finite element model refinement for a structural testbed at NASA Langley Research Center 
(Controls-Structures-Interaction Evolutionary model) is presented to demonstrate the feasibility of the proposed 
approach. 


Introduction: Problem Statement 

T YPICALLY, the spacecraft structure can be modeled as a linear, 
time- invariant flexible system, which in turn can be represented 
by the following second-order dynamical equations: 

Mx + Dx + Kx = Hf (1) 

where M is the positive definite mass matrix, D is the positive 
definite (semidefinite in the presence of rigid-body modes) damp- 
ing matrix, K is the positive definite (semidefinite in the presence 
of rigid-body modes) stiffness matrix, H is the disturbance input 
influence matrix, x is a k x 1 vector of displacements, and / is a 
e x \ vector of disturbances to the system. Usually, a finite clement 
analysis is used to obtain these matrices analytically. However, the 
accuracy of the finite element model in predicting the dynamical 
behavior of the structure depends on a number of factors, such as 
proper knowledge of element and component material and geo- 
metric properties, appropriate meshing, correct joint modeling, etc. 
From past experience with flexible structures, the accuracy of the 
finite element model is limited when compared to test results from 
modal parameter identification. In almost every structure, the modal 
frequencies and amplitudes predicted using finite element models 
differ from those obtained from modal testing. This lack of accuracy 
in modal parameters can be a detriment to control system design. 
Control system design for flexible systems is challenging because 
of their special dynamic characteristics: a large number of struc- 
tural modes within the controller bandwidth; low, closely spaced 
modal frequencies; very small inherent damping; and insufficient 
knowledge of the parameters. 1 

Control system design requires accurate knowledge of the plant 
that is to be controlled. In the case of spacecraft control systems, 
this means that an accurate knowledge of the parameters associated 
with the flexible modes of the spacecraft, such as modal frequencies, 
damping ratios, and mode shapes, is required. The need for accu- 
rate knowledge is particularly critical for the modal frequencies. In 
traditional gain-stabilized spacecraft control design, this knowledge 
is required to achieve nominal performance while guaranteeing sta- 
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bility margins in the form of phase and gain margins. In modem 
control system design, which may be gain or phase stabilized, this 
knowledge is required to achieve nominal performance as well as 
specific degrees ol stability and performance robustness. 

One approach to obtain accurate models of flexible structures is to 
use models that can be extracted directly from system identification 
(ID) or modal test data. This is a teas i ble approach, which has been 
quite successful in a number of applications. However, the useful- 
ness of this approach is limited in that the refined model obtained 
applies only to the hardware configuration of the system ID. In other 
words, the model obtained from system ID data is only valid at the 
input/output channels that are used in the test setup. If the model of 
a component changes, additional inputs or outputs are included, or 
simply new elements are added, the model obtained through system 
ID loses its relevancy unless additional system ID tests are per- 
formed. Moreover, these models do not easily lend themselves for 
other required performance and reliability analyses, such as stress 
and strain analysis, vibration and jitter analysis, etc. To overcome 
the limited aspect ol the system ID models, one can use an analyt- 
ical model obtained through finite element analysis, provided that 
these models can be made to have sufficient accuracy for design and 
analysis. Thus, in this paper we address the problem of refining the 
analytical model of the flexible spacecraft using the system ID data. 

To date, different techniques have been proposed for refining the 
finite clement model of a flexi ble structure based on modal testing or 
system ID procedures. 2 - 3 Model refinement involves techniques that 
refine the finite element model by minimizing the level of disagree- 
ment between the model and test results. These techniques gener- 
ally start with a set of parameters of the model (typically, physical 
parameters at the clement level, e.g., material and geometric prop- 
erties) and systematically tune those parameters to reduce or min- 
imize some measure of disparity between the model and test data. 
This may be a time-related measure, for example, the difference 
in time history responses, or a modal-related measure, for exam- 
ple, difference in modal frequencies. Various optimization schemes 
and least-square techniques have been suggested for the refinement 
process. 

This paper describes a novel approach for the refinement of fi- 
nite element models. The approach presumes that modal analysis 
or system ID tests have been performed and that modal parameters, 
such as frequencies, damping ratios, and mode shapes (at sensor 
locations), have been identified for modes in the range of interest. 
The proposed approach models the possible refinements in the mass, 
damping, and stiffness matrices of the finite element model in the 
form of a constant gain feedback with acceleration, velocity, and 
displacement measurements, respectively. The freedom to change 
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model parameters, as well as the relative degree of change desired 
in one parameter with respect to the rest, is embedded in the el- 
ements of the input and output influence matrices for the various 
measurements. Once the elements of the input and output influ- 
ence matrices have been defined and fixed, the problem of model 
refinement reduces to obtaining position, velocity, and acceleration 
gain matrices that reassign a desired subset of the eigenvalues of the 
model, along with partial mode shapes, from their baseline to target 
values. Hence, the problem of model refinement becomes a problem 
of eigensystem assignment with output feedback. However, symme- 
try and the positive definiteness requirement of the mass matrix, and 
the positive definiteness (semidefiniteness if rigid-body modes are 
present) requirement of the stiffness and damping matrices, neces- 
sitates that gain matrices should be constrained such that the refined 
mass matrix remains symmetric and positive definite and the refined 
stiffness and damping matrices remain symmetric and positive defi- 
nite (semidefinite). Moreover, the refinements in the system matrices 
should accommodate the sparsity of the nominal model, such that 
unwarranted additional connections are not introduced through re- 
finement. In this paper, a procedure for obtaining symmetric gain 
matrices via eigensystem assignment is described first. To perform 
the required eigensystem assignment, a modified procedure to the 
sequential algorithm outlined in Ref. 4 is followed. The modified 
procedure provides the ability to use acceleration feedback, needed 
to refine the mass matrix, as well as the capability to assign eigenvec- 
tors partially. Second, additional constraints, in the form of quadratic 
inequality constraints, are outlined to render the symmetric gain 
matrices dissipative, thus guaranteeing the stability of the refined 
model. Finally, quadratic equality constraints are introduced to ac- 
commodate the sparsity requirements of the refined model. A nu- 
merical example involving model refinement of a structural testbed 
at NASA Langley [CSI evolutionary model (CEM) phase II] is pre- 
sented to demonstrate the application ot this approach. 

Model Refinement 

By observing of the nominal dynamical model of the system given 
in Eq. (1), the dynamics of the refined system may be written as 

(M -F AM)i + (D + A D)x -F (K + A K)x = Hf (2) 

where AM is a symmetric matrix representing the refinement in 
the mass matrix, satisfying (M + AM)>0; AD is a symmetric 
matrix representing the refinement in the damping matrix, satisfy- 
ing (D + AO) > 0; and AK is a symmetric matrix representing the 
refinement in the stillness matrix, satisfying ( K -F AK) > 0. Note 
that the positive definiteness conditions for the refined stillness and 
damping matrices reduce to positive semidefiniteness in the pres- 
ence of rigid-body modes. Now, expand the refinement gain matrices 
as follows: 

AM = L M G M L T M , AD = L D G D L T Di AK = L K G K L T K 

(3) 


where L M is a matrix representing the distribution of refinements 
that are allowed in the mass matrix. The elements of matrix L M can 
vary depending on what elements in the mass matrix are chosen to 
vary. For example, if the chosen element of the mass matrix is the 
one at the / th row and column, then all of the elements of the / th row 
of matrix L M may be chosen to be zeros, except one, which is set to 
1. The matrix Gm represents the symmetric gain matrix associated 
with the mass matrix (acceleration gain matrix) that determines the 
extent of the refinement. The matrices Lp, G D , L K , and G K are 
similarly defined to characterize the refinement for the damping 
and stiffness matrices, respectively. Note that the refinements in the 
mass, damping, and stiffness matrices may be viewed as symmet- 
ric, constant-gain acceleration, velocity, and position feedback. The 
system equations for the refined system may be rewritten as follows 

Mx 4- Dx + Kx = Hf -F « 


u = — \LkGkL\ LpG pL D LmGmL^ 


(4) 


Assume that a number of modes in the desired frequency range 
have been identified via a system ID procedure, and let the modal fre- 
quencies, damping ratios, and modal amplitudes (at the disturbance 
and measurement locations) be denoted by ft n Z {% and 4> f , respec- 
tively. Here, ft, is an r x 1 real vector of natural frequencies of the r 
identified modes; Z, is an r x 1 real vector of modal damping ratios; 
and 4> f is an 5 x r complex matrix, whose r columns represent the 
mode shapes of these identified modes at 5 locations. After noting 
that, for real systems, complex eigenvalues occur in conjugate pairs, 
let the target eigenvalues and eigenvectors be defined as 

a ;‘ ~ 1 = -z;n; + jnifi - zf, - 1 = <*>', 


a? = -zjn; - jti;/ i - zf, *, 2i = <t>; 

1 = 1 r (5) 

The overbar in the expressions in this section refer to complex conju- 
gation of the elements of the corresponding vector (or matrix) only, 
as opposed to the Hermitian operator, which involves transposition 
and complex conjugation. Now the problem of model lefincment 
may be expressed as the problem of finding symmetric accelera- 
tion, velocity, and position gain matrices (G w , Gp, and G K ) such 
that the 2 r eigenvalues of the system 

(M + LmG m L t m )x + (D -F L d G d L d )x 

-F {K 4- LkG k L t k }x = Hf ^ 

are assigned to Aj. / = 1, 2, . . . , 2 r, and the 5 elements of corre- 
sponding eigenvectors are assigned to 4*/ , / = 1 , 2 2r , subject 

to the condition that the refined mass, damping, and stillness ma- 
trices are positive definite. The partial assignment of eigenvectors 
may be defined as 

R<t> = vF, 0) 

where matrix R represents the influence coefficient matrix tor the 
system ID sensor locations. The procedure developed and followed 
to compute the gain matrices is described in the next section. 

Refinements in Damping Matrix 
via Eigensystem Assignment 

The task of assigning the eigenvalues and partial eigenvectors ot 
the system in Eq. (6) with symmetric output feedback gam matrices 
is accomplished using a sequential algorithm described by Ref. 4. 
The algorithm is modified here to accommodate the partial eigen- 
vector assignment and acceleration terms to include refinements in 
the mass matrix. For the simplicity of presentation, the procedure is 
described for refinements in the damping matrix alone and then tor 
an all inclusive model refinement. 

In each step of the sequential procedure, one conjugate pair ot 
eigenvalues is assigned to desired values while making sure that the 
previously assigned eigenvalues are not disturbed. The procedure 
uses a first-order descriptor representation of the system, obtained 
from Eq. (6), 



that is, the descriptor form 

Ez = Az -F Bu + Pf, u = -GpC D z W 

where 



represents the state in the first-order descriptor form, with Cp-B T . 
A brief description of the sequential procedure is given next. The 
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reader is referred to Ref. 4 for a thorough description of the proce- 
dure. 

1) The procedure employs the generalized ordered real Schur 
transformations of the system matrices, E and A. 

2) Orthogonal transformations are used to move previously as- 
signed eigenvalue pairs to the top left block of the pair (E, A), where 
(£, A) are in ordered real Schur form , 5 and the structure of the new 
gain matrix is prescribed such that it only affects the eigenvalues 
in the lower bottom partition of the system matrices. For example, 
assume that (k — 1 ) conjugate pairs of the eigenvalues have been 
placed in the previous steps and that they are in the top left block 
of (E, A). Let Np k denote a matrix whose columns form an orthog- 
onal basis for the left null space of C Qkv that is, N Dk is a matrix 
with orthogonal columns such that A^ C Dk , = 0. Here, C Dkl denotes 
the first 2 (k — l) column partition of the output influence matrix in 
transformed coordinates. If the gain matrix in the transformed co- 
ordinates Gp k is constructed as 

Gp k = Np k G Dk Ny l ( 10 ) 

where G Dk may be an arbitrary matrix, then output feedback with 
the gain matrix will not affect the (k - 1 ) eigenvalue pairs assigned 
in the earlier steps. 

3) At each step, an intermediate gain matrix is computed to as- 
sign a pair of eigenvalues to desired values in lower bottom partition 
of the system matrices. The algorithm used in this eigenvalue as- 
signment was initially developed in Ref. 4. but is modified here to 
accommodate partial eigenvector assignment as well. 

4) At each step, after computing the gain matrix that assigns a pair 
of desired eigenvalues, the intermediate state matrix is transformed 
to a generalized Schur form with all earlier assigned eigenpairs in 
the top left block of the updated system matrix. 

5) The overall gain matrix is constructed by accumulating the 
gains from each step. 

6 ) This process can be continued until up to m eigenvalues have 
been assigned to the desired locations, where m denotes the number 
of inputs or outputs. 

Eigenpair Assignment 

This section describes the approach to select output feedback 
gains to assign one pair of complex-conjugate eigenvalues, while 
ensuring that the gain matrix is symmetric and the partial eigenvec- 
tors are as close as possible to their corresponding target vectors. 
Assume that the klh eigenpair is to be assigned. For notational sim- 
plicity, the system matrices will be denoted as E 2 2 , A 22 , B 2 , and 
C 2 , the output feedback gain matrix will be denoted as G, and the 
desired eigenvalue pair will be denoted (A, A). The problem is to 
select a symmetric matrix G, such that (A, A) is a generalized eigen- 
pair of the closed-loop system matrix, (E n , A 22 — B 2 GC 2 ), and the 
eigenvectors are partially assigned to desired values, as given in 
Eq. (7). 

Let 4> be the eigenvector corresponding to the eigenvalue A. The 
generalized eigenvalue problem becomes (A E 22 — A 22 + B 2 GC 2 )4> = 
0. This expression can be rewritten as 



' 4> ' 


■ 0 ■ 

[XE 21 - A u I B 2 ] 

GC 2 <{> 

= r 

_GC 2 0_ 


It is obvious from Eq. (II), that the vector on the right-hand side 
of the expression must lie in the right null space of T. Let A be a 
matrix whose columns form an orthogonal basis for the null space 
of T, that is, rw = 0. Note that unlike an actual control design 
problem where the number of inputs/outputs are usually fixed, we 
may choose the number inputs (parameters that can be changed in 
the model) large enough as to provide the freedom to assign the 
desired eigenvalues and specified elements of the corresponding 
eigenvectors. Although E 22 , A 22 , and B 2 are real matrices, T and 
A are complex matrices because the eigenvalue A is a complex 
scalar. However, to ensure that the gain matrix is real, the eigenvec- 
tor corresponding to the complex-conjugate eigenvalue is chosen 
to be the complex conjugate of <£, that is, 4> is chosen to be the 
eigenvector corresponding to A. 


Because columns of A span the null space of f, it follows that 


' 0 ‘ 

~ A a = 

’ N { ■ 

_GCi4>_ 


. N '- . 


w-herc a is an arbitrary vector of complex elements and the matrices 
A 1 and A 2 are formed by partitioning A compatibly with and 
GC 2 (j). From Eq. ( 12), 4> = N\Ct and GC 2 <fi= /i 2 a, which leads to 


G Ci Ah a = Aha 


(13) 


The eigenassignment problem is now reduced to selecting a such 
that there exists a symmetric gain matrix G satisfying Eq. (13). With 
4> being the eigenvector corresponding to A, real solutions for the 
gain matrix G can be obtained, and the equations can be written out 
to involve only real arithmetic operations as follows. 

_ For the eigenvalue A, w'ith eigenvector </>, the matrix T — 
[A E 22 — A 22 I /L], and _N is a matrix whose orthogonal columns 
span the null space of f . If the arbitrary coefficient vector is chosen 
to be d, the complex conjugate of a, then it follows that 

G C 2 A s d — Aid (14) 


Equation (13) or Eq. (14) can be rewritten as 


CC2[Re(Ah) 


- Im(Ah)] 


Re(a) 

Im(of) 


= [Rc(Ah) 


— Im(A 2 )| 


Rc(or) 

Im(of) 


(15) 


and 

GCi [Im( A] ) 


RcCAO] 


Re (or) 
lm(a) 


= [Im(Ai) 


Re(Ai)) 


Re(a) 

Im(tt) 


(16) 


where Re( ) denotes real part of the argument and Im() denotes 
imaginary part of the argument. In compact form these equations 
are written as 


G W; p = \\p, G W 2 p = lA/> (17) 

where p = (Re(a r ) Im(a r )] r , = Ci|Re(Aj) -Im(Aj )], Vj = 
[Re(A 2 ) -Im(A 2 )], W 2 = C 2 [Im(A, ) Re(A| )], and V> = 
[lm(A 2 ) Re(A>)]. Note that Eq. (17) is a system of quadratic equa- 
tions in the unknown variables, namely, the elements of the gain 
matrix G and the coefficient vector/?. Furthermore, the elements of 
G should be constrained such that G is symmetric, and the solution 
of the system has to yield an eigenvector, for the whole system, T 
that satisfies the partial eigenvector conditions of Eq. (7), that is. 


RU T = 47 * - 1 


(18) 


where R — [R 0] is the corresponding coefficient in the descrip- 
tor form of the system equations, U is the right unitary matrix in 
the generalized Schur form at the klh step (that keeps the closed- 
loop matrices in real Schur form), and 4 * ~ k 1 is the target partial 
eigenvector for the klh step. Note that, if the condition of Eq. (18) 
is satisfied for one of the eigenvectors of the eigenpair, it would 
identically satisfy its complex conjugate. Here, we assume that the 
set of earlier assigned eigenvalues docs not match the remaining 
eigenvalues of the system, either before or after the eigenpair as- 
signment. This mild assumption ensures that the eigenvectors of the 
earlier assigned eigenvaiues/cigenvectors remain unchanged as ad- 
ditional eigenvalues are assigned. Furthermore, this means that the 
eigenvector condition of the type in Eq. (18) can be imposed one 
mode at a time, and once imposed for an eigenvector, it need not 
be rcimposed again. Now, considering the eigenvalue problem of 
the whole system for the eigenvalue being assigned, one can write 



-A n 

0 


A /: | 2 
A E 22 


- A , 2 + B[ GC 2 

— A 2 2 HK B 2 GC 2 



(19) 
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Note that 



Solving for <p in terms of <j>, one obtains 

ip = — (A£ii — An) *(Ai£i2 — A 12 + BiGCj)^ (20) 
and by substituting for GC 2 <t> and from Eq. (12), one has 
<p = N\ce 

<p — —(XE\\ — Ah) _1 ([X£i 2 “ AnJjVj + B\N 2 )ol = 0a (21) 
or 

T = f a = Sa (22) 

Using Eq. (22) into Eq. (18) and expanding and separating the real 
and imaginary parts, one has 

RUlRe(S) — Im(5)]j^j“jj =Re('I/, 2 *- 1 ) 


eigenvalues may be obtained as follows. Denote Vi = V.p, y 2 = V 2 J) . 
X\ = Wip and x 2 = W2P, and let X - [*, jc 2 ] and Y = [y, y 2 ], then 
Eq. (17) is rewritten asGX = K. LetQbean orthogonal matrix, 
such that 

Q t Y = (29) 

where Y\ is a nonsingular 2x2 matrix (otherwise, the problem is 
solved trivially). The matrix 0 can be obtained by QR factorization 
of Y . Now, define Xu *2 as follows: 

?‘l = Q r x (30) 

X7] 

where X\ is a 2 x 2 matrix, and X 2 is a (m - 2) x 2 matrix. Now 
X\ is nonsingular if X| and *2 are linearly independent (otherwise, 
the problem is trivial). Defining Gu = Pj X ^ 1 , it can be seen that 



Therefore, it follows that the matrix G defined as 


Rt/[Im(S) Re(S)] ') (23) 




(32) 


Recalling the definition of vector p from Eq. (17) and combining 
these equations, one obtains 

Lp = q (24) 

where 

RU Re(S) —RU lm(S)l = jRe(>R, 2 *- ') 1 

RUlm(S) /?t/Re(S)J’ Q [lm(4' 2 *- | )j 

Therefore, a coefficient vector p and a symmetric gain matrix G that 
satisfy Eqs. (17) and (24) have to be found. The condition for the 
existence of a symmetric gain matrix G that satisfies Eq. (17) has 
been established in Ref. 4 and is given as the existence of a vector 
p that satisfies 

p T (vfW 2 - vj W,)p=p T Jp = 0 (26) 

To summarize, the conditions for the placement of an eigenpair of the 
system to desired values, while partially assigning the corresponding 
eigenvectors to target values, reduces to computing a coefficient 
vector p that satisfies the quadratic equation given by Eq. (26) and 
the linear system of equations represented by Eq. (24). 

One possible approach to obtaining a coefficient vector p that 
satisfies Eqs. (26) and (24) would be to first solve for p in Eq. (24), 
to obtain 

p = L^q 4 - N l (3 (27) 

where () + is the pseudoinverse of (), N L is a matrix collecting a 
set of basis vectors for the right null space of matrix L, and /3 is a 
coefficient vector associated with the basis vector, yet to be defined. 
Substituting for p from Eq. (27) into Eq. (26) yields 

f(f3) = /3 t N[ J N l (3 + f3 T N[ J L + q + q T L +T J N l {3 

+ q r L+ T JU ¥ q = 0 (28) 

Standard Newton methods for obtaining the solution of nonlin- 
ear equations may be used to obtain a solution vector /3. Ana- 
lytic gradients of f(J3) are readily available because the gradi- 
ent of any quadratic, f(/3) = /3 r 0 i/3 + 20^/3 H-c, is given by (3/ 
= (0i + Q])@ + 202- The nonlinear problem of Eq. (28) 
is very well behaved because the function is quadratic in f3 and an- 
alytic gradients are linear in (3. 

Once a coefficient vector p (or /3) has been determined, the sym- 
metric gain matrix G that assigns the desired complex-conjugate 


satisfies GX = Y 

Note that if only eigenvalue assignment is required, that is, it is not 
required to perform partial eigenvector assignment, then the solu- 
tion vector p does not have to satisfy Eq. (24), and only the quadratic 
equation given by Eq. (26) needs to be satisfied. The solution to this 
equation can be obtained through standard Newton methods as men- 
tioned earlier. However, because Eq. (26) is a simple quadratic, the 
check for existence and computation of a solution can be achieved 
via examination of the matrix J. A solution vector exists if and only 
if the symmetric part of J is either indefinite or semidefinite. If the 
symmetric part of the matrix has zero eigenvalues, then any cor- 
responding eigenvector is a solution for vector p. If the symmetric 
part of J is indefinite, then any linear combination of eigenvectors 
of the symmetric part, whose corresponding eigenvalues are not all 
of the same sign, qualifies as a solution, if the coefficient of the 
linear combination are chosen such that the quadratic in Eq. (26) 
vanishes. 

Once the gain matrix G is computed, the current refinements in the 
damping matrix, represented by Go k is determined from Eq. (10), 
and the overall refinement is updated as 

G d *r- G d + G o k (33) 

The procedure described thus far determines a symmetric gam ma- 
trix that reassigns a desired subset of the eigenvalues ot the model, 
along with partial mode shapes, from their baseline to target val- 
ues. However, the symmetry of the gain matrix does not necessarily 
guarantee that the refined (combined) model remains stable. Be- 
cause in most situations the flexible system is open-loop stable, any 
refinements to the analytical model should be such to maintain that 
stability. One approach to this could be to use the design freedom in 
the solution vector p and impose constraints on eigenvalues of the 
refined damping matrix. However, this could be cumbersome, par- 
ticularly, when the size of the system is large (thousands or hundreds 
of thousands of degrees of freedom). Another approach could be to 
use the freedom beyond eigensystem assignment [quadratic equali- 
ties in Eqs. (26) or (28)] to reduce some measure of the gain matrix 
G that represents the refinement in the damping matrix. This can 
be accomplished, for example, by imposing inequality constraints 
in the form of 

1 P T v? Wjp | > S ijP T V, T V iP , (=1,2, j — 1,2 (34) 

where 6^ are positive scalars whose values can be used to adjust the 
norm of the gain matrix. 
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Alternatively, one can require that the current gain matrix G Dk 
be dissipative at every sequence. Although the dissipativity require- 
ment can be constraining, it will guarantee that the refined system 
remains stable. In other words, at every sequence, a pair of eigen- 
values are assigned via a symmetric and dissipative gain matrix. 
Reference 4 provides a set of constraints to impose dissipativity of 
the gain matrix in this setting. These constraints are in the form of 
quadratic inequality constraints in the solution vector/? as follows: 

/. (P) = P T { vf W, + i ( v, r w 2 + Vi w,) Jp > 0 
Mp) =P T {v[w, - ±(VfW 2 + V 2 t W,)}p > 0 


where W r includes all rows of W except rows j \ , j 2 j, r From 

linear matrix theory, V { is in the range space of W t . iff all rows of 
W r are orthogonal to the right null space of V ih that is, 

W c Vf = 0 (38) 


The f th row of matrix V is a 1 x 2 vector defined as [V u p V 2 ,p], 
where the subscript / refers to the /th row of the matrix. The right 
null vector of V-, may be written as 


V, 1 


VvP 

-Vup 


(39) 


Mp) =p t [v 2 t w 2 + + v 2 t w,) \ p > o 


Using Eq. (39) in Eq. (38) and expanding W r as [W u p W 2( p ], one 
obtains 


Mp) = p T { v 2 t w 2 - |( w 2 + v 2 t Wt ) } P > o (35) 


W u pV 2 iP - w 2l p K IP = o (40) 


These quadratic constraints go well with the quadratic symmetry 
condition given in Eq. (26), and hence, the appealing computational 
nature of the algorithm is retained. Another possible approach could 
be to require that at each sequence of the eigensystem assignment 
procedure the overall gain matrix Go remains positive semidefinite, 
that is, the gain matrix is dissipative. 

Sparsity Accommodation 

In most applications involving flexible systems, the system matri- 
ces are not fully populated matrices. On the contrary, these matrices 
exhibit significant sparsity. However, the formulation presented thus 
far for model refinement does not consider or incorporate possible 
sparsity in the system. It is essential to have such a capability because 
it may not be practical to produce refined models that have different 
sparsity (bandwidth) than the nominal model. One exception can be 
in applications wherein identification of possible missing elements 
(connectivity) in the nominal model is of concern. 

To accommodate sparsity in the damping refinement, assume that 
the refinement influence matrix L D is defined such that each of its 
columns has only one nonzero element. In other words, no coupling 
in matrix Lp is considered between the degrees of freedom con- 
sidered for refinement. This assumption would isolate any potential 
couplings between the degrees of freedom in the elements of the 
gain matrix Gp. Also, this assumption would not result in any loss 
of generality because any coupling in influence matrix Lp can be 
absorbed into the gain matrix. The procedure to incorporate sparsity 
in the damping matrix refinements follows. 

As mentioned earlier, the placement procedure is sequential 
wherein at each sequence the model is refined such that a pair of 
complex-conjugate eigenvalues (and possibly partial eigenvectors) 
are reassigned to match system ID values. At step k of this proce- 
dure, the local gain matrix G must satisfy the conditions established 
in Eq. (17), such that 


Therefore, the condition for existence of a gain matrix G, with 
vanishing super diagonal elements j \ , j 2 , . . . , j (/ in its /th row, is 
a system of quadratic equalities in the coefficient vector/?. Similar 
conditions can be established for vanishing elements in other rows 
of matrix G. These quadratic conditions, together with the quadratic 
equation given by Eq. (26) for symmetry, the quadratic inequalities 
given by Eq. (35) for dissipativity, and the linear system of equations 
provided by Eq. (24) for partial eigenvector assignment constitute 
the complete conditions for sparse damping matrix refinement. 

The computation of the gain matrix in the presence of sparsity 
conditions requires a different treatment than described earlier [see 
Eq. (32)]. Here, the gain is directly computed from the solution of 
simultaneous equations that impose eigenvalue assignment, sym- 
metry, and sparsity conditions. Once, a feasible coefficient vector 
/?, satisfying the quadratic and linear requirements described earlier 
for sparse refinement, has been obtained. Then, the system of simul- 
taneous equations that yields the elements of the local gain matrix G 
are derived as follows. Denote X = [ W\p Wjp] and Y = ( V { p VS/?]; 
then transposing Eq. (36) gives 

X t G = Y t (41) 

Packing the columns of matrix G to obtain vector g, Eq. (41) may 
be rewritten as 

r x t o oi 



0 

0 0 X T _ 

where vector Y f) denotes the columns of matrix Y T packed in a 
columnwise manner. The condition for symmetry of matrix G may 
be written as 


G[W x p W2p] = [V iP Vxp] (36) 


Sg = 0 


(43) 


Note that, in this equation and the following development, the sub- 
script k, denoting the step number, has been omitted to improve the 
clarity of presentation. Moreover, for this refinement not to effect 
the (k — 1) eigenpairs earlier matched by the algorithm, the overall 
gain matrix at step k, G, must be defined via Eq. (10). The required 
sparsity in the refinement matrix should be incorporated in the el- 
ements of matrix G. Now assume that sparsity requires that super 
diagonal elements j [y j 2l . . . , j ( , of the /th row of G must be zero, 
where q is the number of zero-valued elements. It can easily be 
shown that Eq. (36) is satisfied identically if 

NGN T N[Wtp W?p] = N[V lP Vip] 

leading to 

G[W\p W 2 p] = [V l p V 2 p]^>GW = V (37) 

For a gain matrix G, with vanishing super diagonal elements yj, 
j 2 , ... , j q in its /th row, to exist while satisfying Eq. (37), the /th 
row of matrix V % Vj, must be in the range space of the rows of W r , 


where S is a coefficient matrix with —1, 0, or 1 for its elements. 
Finally, the condition for proper sparsity of matrix G or NGN T 
may be written as 

Pg = 0 (44) 

where P is a coefficient matrix whose elements depend on the el- 
ements of the matrix N . Now the gain matrix G (or gain vector g) 
can be computed by solving the system of simultaneous equations 
given by Eqs. (42^44). Note that typically there would a great deal 
of sparsity in the coefficient matrices of Eqs. (42^14) so that sparse 
solvers would be ideal for obtaining the solution. Alternatively, one 
can solve first forg from Eqs. (43) and (44), that is, 

g = N SP 0 (45) 

where the matrix N$p is a set of basis vectors for the right null space 
of matrix 

in 
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and 0 is the coefficient vector yet to be determined. Substitute for 
g from Eq. (45) into Eq. (42) and solve for 6 to obtain the gain 
solution. This alternative may be the best approach because there 
may be times where a completely feasible solution for p cannot be 
obtained. In such a case, its best to guarantee that the refinement is 
feasible by being symmetric and having proper sparsity [Eqs. (43) 
and (44)], and then attempt to match the system ID data as close as 
possible [Eq. (42)]. 

Before ending this section on incorporation of sparsity, it should 
be mentioned that if degrees of freedom for model refinement can 
be chosen such that the refinement matrices take a block-diagonal 
form, then a much easier procedure, one that is very similar to the 
case with no sparsity, may be followed. For example, assume that 
the gain matrix at step k is supposed to have the form 


matrices, G w , Gp and G*, is prescribed such that it only affects the 
eigenvalues in the lower bottom partition of the system matrices. 
For example, assume that {k - l) conjugate pairs of the eigenvalues 
have been placed in the earlier steps and that they are in the top leu 
block of (E, A). Let N Mk , Np k , and N Kk be matrices whosepolumns 
form an orthogonal basis for the left null space of C Mki , C Dir and 
C Kli respectively. Here, Cm u , Cp k{ > and C Kk] are the first 2 (k - 1 ) 
column partitions of the output influence matrices in transformed 
coordinates. If the gain matrices (in the transformed coordinates) 
are constructed as 

Cm, = N Ml G Mi N t Mi , G Di = N Dl G 0( N T Dk 

G Kt =N Kt G Kt N T Kt (50) 


NGN 7 


G, 0 

0 g 2 


(46) 


As stated earlier, it can easily be shown that the condition for eigen- 
pair assignment [Eq. (36)] is satisfied identically if 

NGN T N[W ] p W 2 p] = /V[V,p Vtf] 


which, along with Eq. (46), becomes 
A 0 1 

J| . \N[W iP W 2 p] = N[V tP V.jy ) 

0 g 2 J 

leading to 

GxNAWxp W2p] = NdV\P ViP\ 

G 2 N 2 {W [P W 2 P] = N 2 [V\P ViP\ (47) 

where matrices N\ and A 2 are partitions of A according to rows 
of matrices G x and G 2 > respectively. These equation are almost 
identical to those obtained for the case with no sparsity [see Eq. ( 1 7)]. 
The main difference here is that there are two gain matrices for which 
the conditions of symmetry anchor dissipativity have to satisfied. 


General Model Refinements 

The approach for a total model refinement, which includes re- 
finements in the mass, damping, and stiffness matrices, parallels 
the one described in the preceding sections for the refinement in 
the damping matrix alone. In each step of the sequential procedure, 
one conjugate pair of eigenvalues is assigned to desired values while 
making sure that the earlier assigned eigenvalues are not disturbed. 
The procedure uses a first-order descriptor representation of the 
system, obtained from Eq. (6) 



Uo — —G D L T D x = — G d [o 
u K — —G k L t k x = -G d [L t k 



that is, the descriptor form 


(48) 


Ez = Az + Bu + Pf, 


u = 


Um 

Ud 


Uk 


(49) 


The same sequential procedure is now used to assign the eigenpairs, 
one at a time. As in the damping case, the structure of the three gain 


with G Mk , Gp k , and G Kk arbitrary matrices, then output feedback 
with the gain matrices will not affect the {k— I) eigenvalue pairs 
assigned in the earlier steps. 

The approach to select output feedback gains is described next, to 
assign one pair of complex-conjugate eigenvalues, while ensuring 
that the gain matrices are symmetric and the partial eigenvectors 
are as close as possible to their corresponding target vectors. As- 
sume that the kih eigenpair is to be assigned. For notational sim- 
plicity, the system matrices will be denoted as E 22 , A 22 , 

Cdm B Ki , aR d C* 2 , the output feedback gain matrices G M; , 
G Dk > and G Kk will be denoted, respectively, as G (l , G vy and Gj, 
and the desired eigenvalue pair will he denoted (A, X). The profi- 
ler^ is to select symmetric matrices G iJy G lM and G,/, such that 
(A, A) is a generalized eigenpair of the closed-loop system ma- 
trix, (E 22 + B Mi G u Cm 2 , A 22 — Bp 2 G v Cp 2 ~ B/c 2 GjC k 2 ), and the 
eigenvectors arc partially assigned to desired values, as given in 
Eq. (7). 

Let <j) be the eigenvector corresponding to the eigenvalue A. The 
generalized eigenvalue problem becomes 

[}.([■: n + B M2 G a C Ml ) - An + «o,G„C c , + fl*,G, / fy,]0 = 0 

( 51 ) 


This expression can be rewritten as 


£A/i22 — A 2 2 A Bm 2 B p 2 


Bk 2 ] 


<p ' 
G u C M2 <t> 
G v c d , 4> 
G d C Kl <f> 


= r 


<t> 

G a c M2 4> 

G v Cp 2 <p 


= 0 


(52) 


LGrfCjr^J 

It is obvious from Eq. (52), that the vector on the right-hand side of 
the preceding expression must lie in the right null space of T. Let 
A be a matrix whose columns form an orthogonal basis for the null 
space of T, that is, T A =0. Because columns of A span the null 
space of T, it follows that 


' <A ' 


-A/f 

G (J Cm 2 $ 

= Ncx = 

n 2 


G v C Dl <j> 


Ny 

_G d C Kl <t>_ 




where a is an arbitrary vector of complex elements and the matri- 
ces Ai, A 2 , A 3 , A 4 are formed by partitioning A compatibly with 
4>, G u Cm 2 4> , G v Cp4> , and GjCk 2 4>. From Eq. (53), one has 

0 = Aia, G {t C Ml <j> = A 2 a 

G v Cp 2 4> = A 2 a, G d C Kl <t> — A 4 q: (54) 

or 

G (1 Cm 2 N\Ol = A 2 a, G v Cp 2 N\oc = Nycx 


GuCk 2 N x ol ~ N 4 ol 


(55) 
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Following algebraic manipulations, similar to those outlined in 
Eqs. (13-17) for the damping refinement case, the eigenpair as- 
signment problem reduces to the solution of 

G a W Mi p = V W| p, G a W M j> = V M j > , G v W D[ p = V Dl p 

G v W DzP = V D2 P , G, U^p = V Kl p t G d W Kl p = V Kl p 

(56) 


where p = [Re r (a) Im r (a)] r . The matrices W Mi , Wm 2 , W D] , 

, VT* 2 , Vjh,, V„ 2 , VV>, , V D ,, Vk { , and V* 2 are formed from 
the imaginary and real parts of the matrices CV,, Co 2 , CV 2 , 
A 2 , A^ 3, and jV 4 , similar to what was done for the damping re- 
finement case in Eqs. (15) and ( 1 6). Note that Eq. (56) is a system of 
quadratic equations in the unknown variables, namely, the elements 
of the gain matrices G a , G v , and G d and the coefficient vector p. 
The elements of the gain matrices should be constrained such that 
they are symmetric, and the solution of the system has to yield an 
eigenvector T for the whole system, which satisfies the partial eigen- 
vector conditions of Eq. (18) for the kih eigensystem assignment. 
The assumption that the set of earlier assigned eigenvalues does not 
match the remaining eigenvalues of the system, either before or af- 
ter the eigenpair assignment still holds here. Now, considering the 
eigenvalue problem of the whole system for the eigenvalue being 
assigned, one can write 


r x n 

Xnl f 

l\ 

[Ol 






I 

, — / 


X 

j 1 == A. E\ 1 

— A\\ 

L° 

X22JI 

4>\ 

|oj 





12 = ^(£|2 

+ 

,G. 

|Ga* 2 ) — 

^12 

+ B Dl 1 

G v Cp 2 -f 

B K 2i G aC k 2 

22 ~ ^(^22 

+ Bm 22 G u 

,c* 2 ) — 

A 22 

+ Bp 22 

G v C Dl + 

B Kv _G a C K2 


(57) 


Note that 



Solving for <p in terms of </>, one obtains 

<p = -X 2l 'X n 4> (58) 


and by using Eq. (54), one has 


(p — N\a 

ip — —X 1/ ([A £12 — A 12 ] -f N 2 + Bp 2i N3 

+ B K2t N,)a = Qa (59) 


or 


T = 


Q 

iV, 


a = Sol 


(60) 


Using Eq. (60) into Eq. (18) and expanding and separating the real 
and imaginary parts, one obtains an expression similar to the one 
for the damping case [see Eq. (24)] 


Lp=q (61) 

where the matrix L and vector q have been defined in Eq. (25). The 
condition for the existence of symmetric gain matrices G ay G v , and 
G d that satisfies Eq. (56) reduces to the existence of a vector p that 
satisfies 


p T {K t mv, - K, = P T hp = o 

P T ( Wd 2 - Vl Wo, )p=p t J2P = 0 

P T K - vl Wk 1 )p=p t J>p = o (62) 


To summarize, the conditions for the placement of an eigenpair of 
the system to desired values, while partially assigning the corre- 
sponding eigenvectors to target values, reduces to computing a co- 
efficient vectorp that satisfies the three quadratic equations given by 
Eq. (62) and the linear system of equations represented by Eq. (6 1 ). 
This is very similar to the problem obtained for the damping refine- 
ment case, with the exception that instead of one quadratic equation 
we have three quadratic equations. Hence, the approach proposed 
for the damping case, which involved a combination of the solu- 
tion of the linear system of equations along with standard Newton 
methods, may be used to solve for a feasible coefficient vector p. 
Once a coefficient vector p is obtained, the procedure to compute 
the gain matrices G a , G,„ and G d is straightforward and follows the 
treatment described for computing the gain matrix in the damping 
refinement case [see Eqs. (27-32)]. Once the gain matrices G a , G v , 
and G d (or matrices G Mk , Gp k ., and G Kk ) are computed, the current 
refinements in the mass, damping, and stiffness matrices, repre- 
sented by G Wi , G p K , and G Kk , are determined from Eq. (50), and 
the overall refinements are updated as 

G m G M -f G Mk , G o G p + G . , G k G k + G K k 

(63) 

The procedure outlined determines symmetric gain matrices G M , 
G Dy and G K that reassigns a desired subset of the eigenvalues of the 
model, along with partial mode shapes, from their baseline to target 
values. As described for the case of damping matrix refinements, 
the symmetry of the gain matrices does not necessarily guarantee 
that the refined (combined) model remains stable. Because, in most 
situations, the flexible system is open-loop stable, any refinements 
to the analytical model should be such to maintain that stability. One 
approach to this could be to use the design freedom in the solution 
vector p and impose constraints on eigenvalues of the refined mass, 
damping, and stiffness matrices. However, this could be cumber- 
some, particularly when the size of the system is large (thousands 
or hundreds of thousands of degrees of freedom). Another approach 
could be to use the freedom beyond eigensystem assignment to re- 
duce some measures of the gam matrices G^, Gp, and G* that 
represent the refinement in the system matrices. This can be accom- 
plished, for example, by imposing inequality constraints similar to 
those presented in Eq. (34) for each of the matrices. 

Alternatively, one can require that the current gain matrices G M k , 
G p K , and G Kk be dissipative at every sequence. Although dissipativ- 
ity requirement can be constraining, it will guarantee that the refined 
system remains stable. In other w-ords, at every sequence, a pair of 
eigenvalues are assigned via a symmetric and dissipative gain ma- 
trices. Similar to the damping case, dissipativity of the gain matrices 
can be achieved through a set of 12 (4 per gain matrix) quadratic 
inequality constraints in the solution vector p The form of the in- 
equality constraints for each gain matrix is exactly the same as the 
ones given in Eq. (35), except that the appropriate coefficient ma- 
trices are used instead of matrices W W 2 , V\ , and V 2 . Yet another 
approach could be to require that at each sequence of the eigen- 
system assignment procedure the overall gain matrices G m , Gp, 
and Gk remain positive semidefinite, that is, the gain matrices arc 
dissipative. 

Finally, the conditions for accommodating sparsity requirements 
of each of the system matrices directly follows the development 
for the damping refinement. These conditions are expressed as 
quadratic equalities in the coefficient vectorp, similar to those given 
in Eqs. (40) or (47). 


Numerical Example 

The approach for model refinement using eigensystem assign- 
ment has been applied to a finite clement model of the phase II 
CEM, a testbed for control of flexible space structures at NASA 
Langley Research Center. Here, the proposed approach is used to 
refine the damping and stiffness matrices of the structure using sim- 
ulated identified modal frequencies and damping ratios. 

The phase II CEM structure consists of a 62-bay central truss 
(each bay is 10-in. long), along with two horizontal booms for 



690 


MAGHAMI 


suspension, a vertical laser, and a vertical reflector tower, as shown 
in Fig. 1 . This structure has 10 modes with frequencies up to about 
5 Hz and 95 modes with frequencies under 60 Hz. The first six modes 
are rigid-body modes, due to suspension of the structure from the 
laboratory ceiling, that have frequencies up to about 0.3 Hz. Eight 
control stations housing collocated and compatible sensors and actu- 
ators are located at the bays shown in Fig. 1 . Air thrusters, providing 
linear forces, are available at these locations along the directions 
shown in Fig. 1. Linear velocities are assumed to be available at 
these locations along the same directions. 

An eight-degree-of-freedom structural model, which includes the 
first eight modes of the structure, is obtained following dynamic 
condensation techniques and is used in this numerical example. A 
low inherent damping ratio of 0.1% has been assumed for the each 
of the eight modes. The nominal eigenvalues along with damping 
and frequencies are shown in Table 1. Assume that only modes 2 
and 8 are to be considered for refinement and that the frequency of 
mode 2 is low by 10% and its damping ratio is low by almost 25%, 
and the frequency of mode 8 is high by 12% and its damping ratio 
is low by almost 10%. Moreover, assume that the mass matrix is 
perfectly known, such that no refinements in the matrix is required. 
However, it is desired to refine the damping and stiffness matri- 
ces, using the proposed eigensystem assignment technique, such 
that the frequencies and damping ratios of modes 2 and 8 of the 
refined system matches the identified values. No identified eigen- 
vectors are included, that is, there is no need for partial eigenvector 
assignment. 

Assume that there is uncertainty in the elements of the damping 
and stiffness matrices corresponding to degrees of freedom 1, 2, 7, 
and 8. The input/output influence matrices [see Eq. (3)] were then 
chosen as 


L k 


10 0 0 
0 10 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 ’ 
0 0 0 0 
0 0 10 
0 0 0 1 


4 = 0.1x L k (64) 


to provide freedom for the appropriate elements of damping and 
stiffness matrices to change. 


Symmetric Refinements 

The first objective was to decrease the natural frequency of the 
second mode by 10% and increase its damping ratio to 25% so that 
the first target pair of eigenvalues was Xj ,2 = —0.228 ± 0.8841 j. 


First, symmetric position and rate gain matrices were sought to re- 
assign this pair of eigenvalues. Following the procedure described 
in the preceding section, the gain matrices were computed from the 
solution of system of quadratic equation given in Eq. (62), except 
for no equations corresponding to mass matrix refinements. The 
quadratic equations were solved using the MATLAB® (Ref. 5) non- 
linear equation solver routine FSOLVE, which uses a Levenberg- 
Marquardt method. The eigenvalues of the system, with the inter- 
mediate position and rate gain matrices in place, are provided in 
Table 2. Table 2 indicates that the complex-conjugate pair were 
successfully reassigned to desired values. However, the resulting 
refined system has an unstable pole on the real axis. This is to be ex- 
pected because, as mentioned earlier, the symmetry of gain matrices 
does not typically guarantee the stability of the system. 

For the second step, the damping ratio in the mode 8 of the system 
was to be increased to 10%, while its frequency was to be decreased 
by 12%, resulting in the second pair of desired eigenvalues to be 
A 3 ,4 = -1.3 149 ± 13.0835;. Following the sequential approach out- 
lined, first the pair of complex-conjugate eigenvalues were placed 
on the top left partition of the real Schur form of the system using 
orthogonal transformations. 6 Then, the gain matrices were defined 


Table 1 Nominal eigenvalues 


Open-loop 

Damping 

Frequency, 

eigenvalues 

ratio 

rad/s 


-0.0008 ±0.81 80; 
-0.0008 ±0.8301 ;' 
—0.0009 ± 0.8565 ; 
-0.0011 ± 1.1308; 
-0.001 1 ±1.1401 ; 
-0.0019± 1.9100; 
-0.0107 ± 10.7278; 
-0.0149 ± 14.9425; 


0.0010 

0.8180 

0.0010 

0.8301 

0.0010 

0.8565 

0.00 10 

1.1308 

0.0010 

1.1401 

0.0010 

1.9100 

0.00 10 

10.7278 

0.0010 

14.9425 


Table 2 Eigenvalues of refined system 
with intermediate gains 


Closed-loop 

eigenvalues 

Damping 

ratio 

Frequency, 

rad/s 

-0.5018 ±0.6423; 

0.6156 

0.8151 

-0.0068 ±0.8404; 

0.0081 

0.8404 

-0.1 148 ±0.8554; 

0.1330 

0.8630 

-0.0068 ±0.8905; 

0.0076 

0.8906 

—0.2283 ±0.8841 j 

0.2500 

0.9131 

-0.0222 ± 10.6262; 

0.0021 

10.6262 

— 0.3956 ± 14.4418; 

0.0274 

14.4472 

6.2235 

- 1 .0000 

6.2235 

-19.3636 

1.0000 

19.3636 





Y 


Fig. 1 Schematic of phase II CEM structure, with location of eight control stations. 
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Table 3 Eigenvalues of refined system 


Closed-loop 

eigenvalues 

Damping 

ratio 

Frequency, 

rad/s 

—0.4955 ± 0.57 15; 

0.6551 

0.7564 

-0.1219 ±0.8600y 

0.1403 

0.8685 

-0.0054 ±0.881 3; 

0.0062 

0.8813 

-0.2283 ±0.8841; 

0.2500 

0.9131 

-0.1458 ±4.9129; 

0.0297 

4.9151 

— 0.1640± 12.8917; 

0.0127 

12.8927 

-1.3149 ± 13.0835; 

0.1000 

13.1494 

-5.3957 ± 15.7620; 

0.3239 

16.6600 


such that the eigenpair remain unchanged [^ee Eq. (50)] while the 
new pair of eigenvalues were assigned. The cumulative position 
and rate gain matrices Go and G Ky which assign the two pairs of 
complex conjugate eigenvalues, were computed to be 

'30.367 78.615 68.806 24.593' 

G = 78.615 369.249 220.096 63.317 

68.806 220.096 233.235 94.726 

_ 24. 593 63.317 94.726 40.439. 

"134.819 -26.937 -27.664 - 81 . 395 " 

r -26.937 15.249 23.763 9.922 

Gk = (65) 

-27.664 23.763 38.679 6.733 

81.395 9.922 6.733 36.126. 

The eigenvalues of the system, with the position and rate gain matri- 
ces in place, are provided in Table 3, where it is observed that the two 
pairs of eigenvalues had been successfully assigned to the identified 
values. Furthermore, the resulting refined system is stable, although 
there were no measures imposed to guarantee such stability. Also, 
note that the remaining eigenvalues (those that were not reassigned) 
have changed, some significantly. One could make some of those 
eigenvalues invariant during refinement by allowing more elements 
of the damping and stiffness matrices to change. The refinements 
in the damping and stiffness matrices are computed from Eq. (3). 
From these equations, the refinements in the damping and stiffness 
matrices, namely, A D and A K , are of the same order as the matrices 
themselves. However, because of the structures of the assumed L D 
and L k> only the elements corresponding to degrees of freedom 1, 
2, 7, and 8 are nonzero, and are given as 

AD, =0.01 x G D , AK) = G* (66) 

where Go and G K are given in Eq. (65). 

Dissipative Refinements 

As mentioned earlier, there are a number of ways of guarantee- 
ing that the refined system remains stable. One of the proposed 
approaches was to take advantage of the freedom beyond eigensys- 
tem assignment and to determine the solution vector p such that the 
gain matrices, representing the refinements in the model, are dissi- 
pative. In the second example, the same model refinement problem 
as in the first case was considered with the exception that position 
and rate gain matrices were constrained to be dissipative. The po- 
sition and rate gain matrices were computed from the solution of 
the system of quadratic equalities, given by Eq. (62), and quadratic 
inequalities, given by Eq. (35), except for no equations correspond- 
ing to mass matrix refinements. The gain matrices were determined 
to reassign the eigenvalues of the second mode to its target val- 
ues at A., i2 = -0.228 ±0.8841 j. The system of quadratic equalities 
and inequalities were posed in the form of a minimax problem and 
was solved using MATLAB’s minimax solver routine MINIMAX. 
The eigenvalues of the system, with the intermediate position and 
rate gain matrices in place, are provided in Table 4. Table 4 indi- 
cates that the complex-conjugate pair were successfully reassigned 
to desired values. The remaining eigenvalues were all stable, that 
is, the resulting refined system was stable. This is to be expected 
because the dissipative nature of the gain matrices guarantees the 


Table 4 Eigenvalues of refined system 
with intermediate gains 


Closed-loop 

eigenvalues 

Damping 

ratio 

Frequency, 

rad/s 

-0.0025 ±0.8346; 

0.0029 

0.8346 

-0.0137 ±0.8636; 

0.0158 

0.8637 

-0.2283 ±0.8841 j 

0.2500 

0.9131 

—0.2313 ± 1.1556; 

0.2002 

1.1322 

-0.0127 ± 1.9112; 

0.0066 

1.9113 

-2.3693 ±4.8741 j 

0.4372 

5.4195 

-0.0633 ± 10.8268; 

0.0058 

10.8270 

-0.3471 ± 14.8576; 

0.0234 

14.8616 


Table 5 Eigenvalues of refined system 


Closed-loop 

eigenvalues 

Damping 

ratio 

Frequency, 

rad/s 

-0.0023 ±0.8454; 

0.0027 

0.8454 

-0.0094 ±0.8707; 

0.0108 

0.8707 

-0.2283 ± 0.8841 j 

0.2500 

0.9131 

-0.3734 ± 1 .1556 j 

0.2943 

1.2686 

—0.1416 ±6.3839 j 

0.0222 

6.3855 

- 1.3149 ± 13.0835 ; 

0.1000 

13.1494 

-0.9821 ± 13.2494; 

0.0739 

13.2857 

-0.7522 ± 14.4382; 

0.0520 

14.4577 


stability of the system. Next, the second pair of eigenvalues was 
reassigned to A 3>4 = — 1 .3 149 ± 13.0835 j . Following the sequential 
approach outlined, first the pair of complex-conjugate eigenvalues 
were placed on the top left partition of the real Schur form of the 
system using orthogonal transformations. 6 Then, the gain matrices 
were defined such that the first eigenpair remain unchanged [sec 
Eq. (50)] while the new pair of eigenvalues were assigned. The cu- 
mulative position and rate gain matrices G D and G K , which assign 
the two pairs of complex-conjugate eigenvalues, were computed to 
be 


~ 77.920 
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22.419 

12.532 

24.411 

22.419 
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6.453 
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60.661 
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42.725 

-3.516 

-7.936 

3.756 

-3.516 
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-1.421 

-121.149 

-7.936 

-1.421 

455.059 


The eigenvalues of the system, with the position and rate gain ma- 
trices in place, are provided in Table 5, where it is observed that 
the two pairs of eigenvalues had been successfully assigned to the 
identified values. Furthermore, the resulting refined system is stable, 
as expected. Also, note that the remaining eigenvalues (those that 
were not reassigned) have changed, some significantly. Again, one 
could make some of those eigenvalues invariant during refinement 
by allowing more elements of the damping and stiffness matrices 
to change. The refinements in the damping and stiffness matrices 
are computed from Eq. (3) and are given in Eq. (66), with the gain 
matrices from Eq. (67). 

Comparisons of the refinements in each example indicate that no 
conclusions can be made in regards to the direction or magnitude 
of the computed refinements. In these examples, the computed re- 
finements in damping matrix for the second example are typically 
lower than those obtained for the first example. However, the situa- 
tion is reversed for the refinements in the stiffness matrix. This may 
be attributed to the variability in the solutions of the minimax opti- 
mization algorithms as well as the nonlinear equation solvers, in the 
sense that they may converge to different solutions depending on the 
starting points. In these examples, the starting estimate for the solu- 
tion vector p was randomly chosen, in each example. Conceivably, 
one could attempt to exploit the freedom beyond eigensystem as- 
signment to minimize, in some sense, the norm of the gain matrices 
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to minimize the effective refinement needed for partial model 
matching. 

Conclusions 

This paper presented a novel approach for the refinement of the 
dynamic model of flexible structures using an eigensystem assign- 
ment technique. The approach presumes that modal parameters, 
such as frequencies, damping ratios, and mode shapes (at sensor 
locations), have been identified for modes in the range of interest. 
The proposed approach models the possible refinements in the mass, 
damping, and stiffness matrices of the finite element model in the 
form of a constant gain feedback with acceleration, velocity, and 
displacement measurements, respectively. The freedom to change 
model parameters, as well as the relative degree of change desired in 
one parameter with respect to the rest, is embedded in the elements 
of the input and output influence matrices for the various measure- 
ments. Once the elements of the input and output influence matrices 
have been defined and fixed, the problem of model refinement re- 
duces to obtaining position, velocity, and acceleration gain matrices, 
with appropriate symmetry (and/or dissipativity) and sparsity that 
reassign a desired subset of the eigenvalues of the model, along 
with partial mode shapes, from their baseline values to those ob- 
tained from system identification test data. Hence, the problem of 
mode refinement becomes a problem of eigensystem assignment 
with output feedback. The proposed procedure assigns one conju- 
gate pair of eigenvalues at each step using symmetric (or symmetric 
and dissipative) output feedback gain matrices, while ensuring that 
the eigenvalues assigned in the earlier steps are not disturbed. More- 
over, the procedure provides that original sparsity of the system ma- 
trices in nominal model is preserved. The advantages of the proposed 
approach are that 1) it provides a systematic and computationally 
tractable means for exact model refinement, that is, the refined model 


would match the identified values exactly, without dependence on 
a nonlinear optimizer, and 2) it characterizes the freedom beyond 
model refinement for possible exploitation, which is inherent in the 
elements of the input and output matrices, as well as the elements 
of the position, velocity, and acceleration gain matrices. This free- 
dom may be exploited to minimize the sensitivity of the refined 
model, to minimize global or local changes to the system matrices, 
etc. A numerical example, involving finite element model refine 
ment for a structural testbed at NASA Langley Research Center 
(CEM) was presented to demonstrate the feasibility of the proposed 
approach. 
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